Matched predictions for Higgs production via heavy-quark loops in the SM and beyond 
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The main Higgs production channel at hadron colliders is gluon fusion via heavy-quark loops. 
We present the results of a fully exclusive simulation of gluon fusion Higgs production based on the 
matrix elements for h + 0, 1,2 partons including full heavy-quark loop dependence, matched to a 
parton shower. We consider a Higgs with standard model couplings as well as models where the 
Higgs has enhanced couplings to bottom quarks (6-philic). We study the most relevant kinematic 
distributions, such as jet and Higgs pr spectra and find that matched samples provide an accurate 
description of the final state. For the SM Higgs, we confirm the excellent accuracy of the large heavy- 
quark-mass approximation also in differential distributions over all phase space, with significant 
effects arising only at large pr- For a &-philic Higgs however, the loops have a dramatic impact on 
the kinematics of the Higgs as well as of the jets and need to be accounted for exactly to achieve 
reliable event simulations. 

PACS numbers: 12.38.Cy, 12.38.-t, 13.85. Qk, 14.80.Bn, 14.80.Ec 
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I. INTRODUCTION 

The CERN Large Hadron Collider (LHC) is running 
at a center of mass energy of 7 TeV and it has already 
accumulated several inverse femtobarns of integrated lu- 
minosity per experiment. One of its main goals is to ex- 
plore the details of electroweak symmetry breaking, and 
in particular to establish the existence of a Higgs sector 
of or beyond the Standard Model (SM). 

At the LHC, Higgs boson production mainly proceeds 
via a quantum effect, gluon fusion (GF) This is 

induced by heavy-quark loops, in particular the bot- 
tom and the top quarks, the latter being by far the 
dominant one in the SM. For a not too heavy Higgs 
boson (to/j < 2m t ), and in appropriate kinematic re- 
gions (j>\ < mt), the top quark can be integrated 
out, resulting, to a very good approximation, in a sim- 
ple, non-renormalizable effective field theory, £heft = 
-\^F^ a F^ u (HEFT) v being the Higgs field 

vacuum expectation value and F^ v,a the QCD field ten- 
sor. The next-to-leading order (NLO) QCD correc- 
tions f5r§] were calculated decades ago in both HEFT 
as well in the full SM and found to be very large 
(ct nlo /ct lo ~ 2). This motivated the formidable en- 
deavour of the next-to-next-to-leading order (NNLO) 
QCD calculations, which have been fully evaluated in 
HEFT @-|Tlj]. The exact NNLO calculation involves 
three loop massive diagrams, and is currently out of 
reach. However, recently, the finite top-quark mass ef- 
fects to the total NNLO prediction have been estimated 
through a power expansion and found to have a 

negligible impact on total rates. Soft gluon resummation 
effects have also been studied in HEFT at NNLL [l6l [l7 | . 
On the other hand, it is known that in the hard tails 
of differential distributions or even in special kinemat- 



ics regimes, such as at small-x [18[, loop effects need to 
be accounted for exactly. So far, the recommended best 
predictions for Higgs GF inclusive production rates in 
the standard model [19] are based on the NNLO+NNLL 
results in HEFT, while keeping the heavy quark mass 
dependence at NLO+NLL (20j . 

In Beyond the SM (BSM) theories, GF becomes sen- 
sitive to all colored states in the spectrum with signifi- 
cant couplings to the Higgs (es). Even though the book- 
keeping becomes more involved, as long as such states 
are heavy, an effective field approximation can still be 
used and QCD corrections can be computed as in the 
SM. The only genuine complication arises, not from ad- 
ditional heavy BSM particles, but from the possibility of 
bottom quarks to have enhanced couplings to the scalar 
(or pseudo-scalar) states of the theory. In SUSY, and 
more generally in type II two-Higgs-doublet-model sce- 
narios, this corresponds to a large tan /? scenario (where 
/3 = V1/V2, uj.,2 being the vacuum expectation values 
of the Higgs doublets coupling to down- and up-type 
fermions, respectively). In this case, the HEFT approxi- 
mation cannot be employed and the accuracy of the best 
available predictions goes down to NLO fl9j |. 

Being of primary importance, total rates and Higgs 
kinematic distributions are now quite well predicted and 
also available via public codes such as ResBos [2l[ and 
HqT [H [23j]. Differential p\ distributions accurate to 
NLO can be also obtained via HIGLU [24| as well as via 
HPro (25|, which both keep the exact bottom- and top- 
quarks mass loop dependence and therefore can be used 
also for predictions of scalar Higgs in BSM. However, in 
experimental analyses, it is also crucial to get as precise 
predictions as possible for exclusive observables that in- 
volve extra jets, such as the jet pt spectra and the jet 
rates, at both parton and hadron level. To optimize the 
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search strategies and in particular to curb the very large 
backgrounds, current analyses both at Tevatron and at 
the LHC select 0-,l- and 2-jet events and perform inde- 
pendent analyses on each sample [26j |. The final system- 
atic uncertainties are effected by both the theoretical and 
experimental ones of such a jet-bin based separation, see, 
e.g., Ref. [27|- In the HEFT, fully exclusive parton- and 
hadron-level calculations can be performed by Parton 
Shower (PS) pro grams such as Pythia [11], Herwig [29| 
and Sherpa [30] in the soft and collinear approximation, 
or with NLO QCD codes matched with parton showers: 
via the MC@NLO [U and POWHEG j^HHI methods. 
However, beyond the HEFT, no fully exclusive predic- 
tion has been available so far. The reason is that one 
needs to compromise between the validity of HEFT and 
the complexity of higher loop calculations. It is how- 
ever possible to get full exclusive control at hadron level 
on the complex event topology at the LHC, while still 
reaching approximately NLL accuracy, with the help of 
recent sophisticated matchin g m ethods between matrix 
elements and parton showers [36l 

In PS programs, QCD radiation is generated in the 
collinear and soft approximation, using Markov chain 
techniques based on Sudakov form factors. Hard and 
widely separated jets are thus poorly described in this ap- 
proach. On the other hand, tree-level fixed order ampli- 
tudes can provide reliable predictions in the hard region, 
while failing in the collinear and soft limits. To combine 
both descriptions and avoid double counting or gaps be- 
tween samples with different multiplicity, an appropriate 
matching method is required. Several algorithms have 
been proposed over the years: the CKKW method, based 
on a shower veto and therefore on event re- weighting (36j 
and MLM schemes, based on event rejection [371 l38j|. 

In this work, we report on the first matched simulation 
of Higgs production in gluon fusion that retains the full 
kinematic dependence on the heavy-quark loops, in the 
SM as well as in generic scenarios with enhanced Higgs 
couplings with bottom quarks, which we dub "6-philic 
Higgs" . 

The paper is organized as follows. We first describe 
our methodology. Then we present our results for a SM 
Higgs. We show that the matching procedure provides 
reliable results both at the Tevatron and especially at the 
LHC and that the effects from massive quark loops are 
indeed mild over all phenomenologically relevant phase 
space. The 6-philic Higgs is considered in the following 
section, where it is shown that loop effects must be in- 
cluded exactly. Gluon fusion production is also compared 
to a matrix-element matched sample for bb h, which 
is the dominant production mode in this scenario. We 
draw our conclusions in the last section. 



II. METHOD 

Our study is based on the /c^-MLM and shower- 
kx matching schemes implemented in Mad- 



Graph/MadEvent [H], interfaced with Pythia 6.4 
for parton shower and hadronization. As explained 
below, we find it convenient to include the effects of 
the heavy-quark loop by simply reweighting the events 
generated via tree-level HEFT amplitudes. 

In the fcr-jet MLM matching schemes [13, IHJ, matrix 
element multi-parton events are produced with a mini- 
mum separation kr cutoff of Q^in- F° r every event, the 
final-state partons are clustered according to the kr al- 
gorithm, and the fcr value for each clustering vertex cor- 
responding to a QCD emission is used as renormalization 
scale for as in that vertex. For the central hard 2 — > 1 or 
2^2 process, the transverse mass m\ — -p\ + m 2 of the 
particle(s) produced in the central process is used as fac- 
torization and renormalization scale. Subsequently, this 
event is passed to the Pythia parton-shower generator. 
There, one of two schemes is employed. Either, the final 
partons (after parton showering) are clustered into jets, 
using the kr algorithm with a jet cutoff of Q^ in > Q^lh-i- 
The jets are considered to be matched to the original 
partons if A>r (parton, jet) is smaller than the cutoff Q J n ^* n . 
If any parton is not matched to a jet, the event is dis- 
carded. For events with parton multiplicity smaller than 
the highest multiplicity, the number of jets must be equal 
to the number of partons. We call this scheme the the 
fcr-MLM scheme. Alternatively, no matching between 
shower jets and partons is done. Instead, an event is re- 
tained provided that the hardest emission in the Pythia 
parton shower is below the scale Q^in (or, for events from 
the highest multiplity, below the scale Q^in ° n °f * ne s °ft~ 
est parton in the event). This is called the "shower- kr" 
scheme, and allows for the matching scale QjL n to be 
set equal to the matrix element cutoff scale Q^in- The 
two matching schemes have been shown to give equiva- 
lent results [40], but for the case of 6-philic Higgs, the 
shower-A-r scheme allows for lower matching scales and 
it is therefore more efficient. 

In order to take into account the full kinematic depen- 
dence of the heavy quark loop in Higgs production, the 
full one-loop amplitudes for all possible subprocesses con- 
tributing to h + 0, 1, 2 partons have been calculated. An- 
alytic expressions have been generated with FeynArts 3.5 
plj . and manipulated with FormCalc 5.3 [42j . The ten- 
sor integrals have been evaluated with the help of the 
LoopTools-2.5 package [Hj], which employs the reduc- 
tion method introduced in Ref. [43| for pentagons, and 
Passarino-Veltman reduction for the lower point tensors. 
The resulting regular scalar integrals are evaluated with 
the FF package (44|. We have also implemented the 
reduction method for pentagon tensor integrals as pro- 
posed in Ref. [45[ for better numerical stability. The 
codes have been used and validated against known re- 
sults in a previous study [46j . The final implementa- 
tion of the calculation includes the contributions from 
bottom and top quarks and their (destructive) interfer- 
ence. The evaluation of multi-parton loop amplitudes 
is, in general, computationally quite expensive. More- 
over, in the case of inclusive matched samples, an ef- 



3 





Cross section 
Higgs mass [GeV] Tevatron LHC @ 7 TeV 


gg -> ft (SM) 
gg -> ft (SM) 
gg — >• ft (6-only) 
bb^h 


140 0.672 pb 12.2 pb 
500 0.003 pb 0.869 pb 
140 3.0 fb 56 fb 
140 4.55 fb 135 fb 



TABLE I: Reference values for total cross sections for Higgs 
production in the SM and considering only 6-loops, used for 
the normalization of the inclusive samples. Results have been 
obtained via the HNNLO [13] and bbh@NNLO [H| codes, 
with mt = 173 GeV, mt = 4.6 GeV, /j,r — fiF = and 
employing the MSTW2008NNLO pdf set [1|. 



ficient event generation needs as a first phase a rather 
thorough exploration of the phase space. It therefore be- 
comes advantageous to devise a method where the map- 
ping of the integrand can be done in a quick (though 
approximate) way and the evaluation of loops limited to 
a small number of points. Our strategy is as follows. 
Parton level events for h + 0, 1, 2 partons are generated 
via MadGraph/MadEvent in the HEFT model, with 
scale choices optimized for the subsequent matching pro- 
cedure. Before passing them to the PS program, events 
are reweigthed by the ratio of full one-loop amplitudes 
over the HEFT ones, r = |7W LO op| 2 /|A^heft| 2 - The 
reweighted parton-level events are unweighted, passed 
through Pythia and matched using the fc^-MLM or the 
shower-Zc-r scheme. All steps are automatic. To vali- 
date the matching procedure, the effect of changing the 
matching cutoff and other parameters such as Q^ ln and 
Qj^P on several distributions, including the n — » n — 1 
differential jet rates have been extensively assessed. 

Finally, we recall that even though matrix elements for 
up to two final states partons are included in the simula- 
tion, the accuracy of the overall normalization of the in- 
clusive sample is only leading order, exactly as in a purely 
parton-shower result. It is therefore legitimate and con- 
sistent to adjust the overall normalization to the best 
available fully inclusive prediction for the corresponding 
process. To this aim, NNLO cross sections (in fact, just 
NLO for a 6-philic Higgs) at the Tevatron and the LHC 
for the scenarios described below have been obtained via 
publicly available codes and collected in Table fl] 



III. 



SM HIGGS PRODUCTION 



To illustrate the results of our simulations for the Teva- 
tron and the LHC at 7 TeV for a standard model Higgs 
we show a few relevant observables in Figs.[TJ(3] We define 
jets via the kr algorithm, with the distance measure be- 
tween parton i and beam B, or partons i and j as k T B = 
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FIG. 1: SM Higgs pr distributions for m h — 140 GeV in gluon 
fusion production at the Tevatron. Results in the HEFT and 
with full loop dependence (LOOP) are compared to the de- 
fault Pythia implementation (which accounts for 2 — > 2 ma- 
trix element corrections) and to the NNLO+NNLL results as 
obtained by HqT [22l [2a|. Curves normalized to the corre- 
sponding total cross sections of Table [I] 



p T , kl^ = min (pt'Pt) \/2(coshAyjj — cos Afaj) / D . 
Here y is the rapidity and <fi is the azimuthal angle around 
the beam direction. The resolution parameter is set to 
D = 1. Jets are required to satisfy \rij\ < 4.5 and 
p T > 30 GeV. For sake of simplicity, we adopt Yukawa 
couplings corresponding to the pole masses, i.e., for the 
top quark m t = 173 GeV and for the bottom-quark 
mass m,b = 4.6 GeV. Other quark masses are neglected. 
Throughout our calculation, we adopt the CTEQ6L1 par- 
ton distribution functions (PDFs) [50] with the core pro- 
cess renormalization and factorization scales ji r = [if = 

m T = <J (pj 1 ) 2 + to^. For the matching performed in 
MadGraph/MadEvent, the fcr-MLM scheme is cho- 
sen, with Q^P = 30 GeV and Q^ m = 50 GeV. 

In Figs. Q] we show the Higgs pt distribution for Stan- 
dard Model Higgs GF production at the Tevatron with 
nth = 140 in a range of pr relevant for experimental anal- 
ysis. We compare matched results in the HEFT theory 
and in the full theory (LOOP) with Pythia with 2^2 
matrix element corrections. We also include the predic- 
tions from the analytic computation at NNLO+NNLL 
as obtained by HqT [22|, [23j|. The curves are all nor- 
malized to the NNLO+NNLL predictions. The three 
Monte-Carlo based predictions agree very well in all the 
shown range of pt, suggesting that for this observable, 
higher multiplicity matrix element corrections (starting 
from 2 — > 3) and loop effects are not important. This 
is the case also for jet px distributions (not shown) in 
the same kinematical range. The NNLO+NNLL predic- 
tion, on the other hand, suggests a slightly softer Higgs 
spectrum. 

In Figs. [HO we show the Higgs and jet px distribu- 
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FIG. 2: SM Higgs p T distributions for m h = 140 GeV and 
rrih = 500 GeV in gluon fusion production at 7 TeV LHC. 
In the upper plot results in the HEFT and with full loop 
dependence (LOOP) are compared over a large range of pr 
values to the default Pythia implementation, which accounts 
for 2 — s> 2 matrix element corrections. In the lower plot the 
low-pr range is compared to the NNLO+NNLL results as 
obtained by HqT [22, 23]. Curves normalized to the corre- 
sponding total cross sections of Table [U 



tions for Standard Model Higgs GF production at the 
7TeV LHC with m h = 140 and 500 GeV. Once again the 
Monte-Carlo based results agree well with each other. As 
expected, loop effects show a softening of the Higgs pr, 
but only at quite high pr- We also see that the heavier 
the Higgs, the more important are the loop effects. This 
is expected, since the heavy Higgs boson can probe the 
internal structure of the top-quark loop already at small 
Pt- The jet pr distributions do confirm the overall pic- 
ture and again indicate loop effects to become relevant 
only for rather high values of the pr- 

The agreement, on the other hand, of the 
NNLO+NNLL predictions at small pr for both Higgs 
masses it is quite remarkable. In this respect, our anal- 
ysis strongly motivates the use of matched samples for 
simulating GF Higgs production at the LHC. Key distri- 



FIG. 3: Jet pr distributions for associated jets in gluon fu- 
sion production of mj, = 140 GeV and mu = 500 GeV Higgs 
bosons at 7 TeV LHC. 



butions, such as the pr of the Higgs, do agree remark- 
ably well with the best available predictions, for example 
NNLO+NNLL at small Higgs pr, and offer improved and 
easy-to-use predictions for other key observables such as 
the jet rates and distributions. In addition, for heavy 
Higgs masses and/or large pr, loop effects, even though 
marginal for phenomenology, can also be taken into ac- 
count in the same approach, if needed. 



IV. B-PHILIC HIGGS PRODUCTION 

In this section we present the results of a simula- 
tion of a 6-philic Higgs. Parameters are the same as 
in the previous section, except that, as explained be- 
low, the top Yukawa coupling is set to zero and the 
matrix-element matching in MadGraph/MadEvent is 
performed through the shower-fcT matching scheme with 

<9match = 10GeV. 

In Fig. |U we show the p\ distributions for GF pro- 
duction at the 7 TeV LHC of a 6-philic Higgs with 
rrih = 140 GeV. We remind the reader that in our cal- 
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FIG. 4: fo-philic Higgs pr distribution a the Tevatron and the 
LHC with m h = 140 GeV. Results in the HEFT approxima- 
tion (red curve) and with full loop dependence (green) are 
shown. Spectrum of Higgs produced via bb fusion in the five 
flavor scheme is also shown. All samples are matrix-element 
matched with up to two partons in the final state. Curves 
normalized to the corresponding total cross sections of Ta- 
ble HI 



culation the bottom-quark and top-quark masses can be 
chosen independently as well as the value of the cor- 
responding Yukawa couplings. We can therefore study 
the production of scalars with arbitrary couplings to the 
heavy quarks such as those appearing in a generic two 
Higgs doublet model or in the minimal supersymmctric 
standard model. For the sake of illustration, we define a 
simplified scenario where the Higgs coupling to the top 
quark is set to zero. In so doing, we study the Higgs and 
jet distributions relative to a "large tan /3" scenario with 
bottom-quark loops dominating. Note that for simplic- 
ity we keep the same normalization as in the standard 
model, i.e., yb/V^ — nib/v with m& — 4.6 GeV, as the 
corresponding cross sections in enhanced scenarios can 
be easily obtained by rescaling. 

In the 6-philic Higgs production, the particle running 
in the loop is nearly massless, and there is no region in 
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FIG. 5: Jet rates for 6-philic Higgs production (where only 
the coupling to the bottom quark is included) . m,h = 140 GeV 
at the Tevatron and LHC (-/i = 7 TeV). 



mh or pt where an effective description is valid. This 
also means that a parton-shower generator alone has no 
possibility of correctly describing the effects of jet radi- 
ation, and genuine loop matrix-elements plus a matched 
description are needed for achieving reliable simulations. 

In fact, the largest production cross section for a b- 
philic Higgs does not come from loop induced gluon fu- 
sion, but from tree-level bb fusion. Phenomenologically, 
it is therefore very important to be able to also generate 
events for this kind of process, which typically leads to 
final states with more 6-jets than the GF production. We 
do so by matching tree-level matrix elements for h+0, 1, 2 
partons (with a hbb vertex) in the five flavor scheme to 
the parton shower. In so doing we provide a complete and 
consistent event simulation of inclusive Higgs production 
in a 6-philic (or large tan 0) scenario. We note in passing 
that a four- flavor scheme, i.e., starting from the leading 
order process gg — > bbh, could also be employed. While 
this latter approach has some important advantages, it 
also offers complications with respect to the simpler five 
flavor scheme. A detailed comparison between the two 
approaches, which are known to be compatible at the 
level of total cross sections (see, e.g., [l9| and references 
therein), would certainly be welcome. Being beyond the 
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scope of this paper, however, we leave it to future work. 

Fig. [5] shows jet rates for the Tevatron and 7 TeV LHC 
for a 6-philic Higgs for two minimum jet p T scales, 30 and 
100 GeV. As is readily seen from the figure, the effect of 
properly including loop effects is significant already with 
a jet p T cutoff at 30 GeV, and increasingly important for 
larger cutoff values. This immediately translates to the 
effect of a jet veto with a given p J T cutoff for the veto. 

V. CONCLUSIONS 

In summary, we have presented the first fully ex- 
clusive simulation of gluon fusion inclusive Higgs pro- 
duction based on the exact one-loop matrix elements 
for h + 0, 1, 2 partons, matched to Pythia parton 
showers using multiple matching schemes implemented 
in MadGraph/MadEvent. We have compared the 
loop reweighted matched results with the corresponding 
HEFT results, Pythia results, and, when possible, with 
NNLO+NNLL predictions. We have considered both the 
SM Higgs and the case of scalar particles with enhanced 
couplings to bottom quarks and studied the most rele- 
vant kinematic distributions, such as jet and Higgs px 
spectra. Our results highlight the relevance of a com- 
plete loop calculation at large px for a standard model 
Higgs and in all phase space for 6-philic Higgs. Such 



improved simulations might be particularly relevant in 
searches performed via multivariate analysis techniques 
where details about the kinematic distributions of the 
Higgs decay products and accompanying jets can have 
significant impact on the results. 

We conclude by stressing that the method employed in 
this work, i.e., using tree- level amplitudes based on an 
effective theory to generate parton-level events and then 
reweighting them by the exact loop amplitudes before 
matching to the shower, is completely general and can 
therefore be applied to any loop-induced process. Work 
towards the automatization of this approach in Mad- 
Graph 5 [5l| via MadLoop [HJ is in progress. 
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